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Abstract:  Traditionally,  robot  control  research  has  focused  on  the  position  tracking  problem 
where  the  objective  is  to  force  the  robot’s  end-effector  to  follow  an  a  priori  known  desired  time 
dependent  trajectory.  Motivated  by  task  objectives  that  are  more  effectively  described  by  on-line, 
state-dependent  trajectories,  two  adaptive  tracking  controllers  are  developed  in  this  paper  that 
accommodate  on-line  path  planning  objectives.  An  example  adaptive  controller  is  first  modified  to 
achieve  velocity  field  tracking  in  the  presence  of  parametric  uncertainty  in  the  robot  dynamics.  The 
development  aims  to  relax  the  typical  assumption  that  the  integral  of  the  velocity  field  is  bounded 
by  incorporating  a  norm  squared  gradient  term  in  the  control  design  so  that  the  boundedness 
of  all  signals  can  be  proven.  An  extension  is  then  provided  that  targets  the  trajectory  planning 
problem  where  the  task  objective  can  be  described  as  the  desire  to  move  to  a  goal  configuration 
while  avoiding  known  obstacles.  Specifically,  an  adaptive  navigation  function  based  controller  is 
designed  to  provide  a  path  from  an  initial  condition  inside  the  free  configuration  space  of  the  robot 
manipulator  to  the  goal  configuration.  Experimental  results  for  each  controller  are  provided  to 
illustrate  proof  of  validation  of  the  approaches. 

1  Introduction 

Traditionally,  robot  control  researchers  have  focused  on  the  position  tracking  problem  where  the 
objective  is  to  force  the  robot  to  follow  a  desired  time  dependent  trajectory.  Since  the  objective  is 
encoded  in  terms  of  a  time  dependent  trajectory,  the  robot  may  be  forced  to  follow  an  unknown 
course  to  catch  up  with  the  desired  trajectory  in  the  presence  of  a  large  initial  error.  For  example, 
several  researchers  have  reported  the  so  called  radial  reduction  phenomena  (e.g.,  [20],  [21])  in 
which  the  actual  path  followed  has  a  smaller  radius  than  the  specified  trajectory.  In  light  of  this 
phenomena,  the  control  objective  for  many  robotic  tasks  are  more  appropriately  encoded  as  a 
contour  following  problem  in  which  the  objective  is  to  force  the  robot  to  follow  a  state-dependent 
function  that  describes  the  contour.  One  example  of  a  control  strategy  aimed  at  the  contour 
following  problem  is  velocity  field  control  (VFC)  where  the  desired  contour  is  describe  by  a  velocity 
tangent  vector  [22],  The  advantages  of  the  VFC  approach  can  be  summarized  as  follows.1 

*This  research  was  supported  in  part  by  U.S.  NSF  Grant  DMI-9457967,  ONR  Grant  N00014-99-1-0589,  a  DOC 
Grant,  an  ARO  Automotive  Center  Grant,  and  in  part  by  the  Defence  Advanced  Research  Project  Agency  (DARPA) 
through  the  Space  Naval  Warfare  System  Center,  San  Diego,  contract  N66001-03-R-8043. 

1  See  [4] ,  [20] ,  and  [22]  for  a  more  thorough  discussion  of  the  advantages  and  differences  of  VFC  with  respect  to 
traditional  trajectory  tracking  control. 
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The  velocity  field  error  more  effectively  penalizes  the  robot  for  leaving  the  desired  contour. 
The  control  task  can  be  specified  invariant  of  the  task  execution  speed. 

Task  coordination  and  synchronization  is  more  explicit  for  contour  following. 


The  ability  for  a  velocity  field  to  encode  certain  contour  following  tasks  has  recently  prompted 
researchers  to  investigate  VFC  for  various  applications.  For  example,  Li  and  Horowitz  utilized  a 
passive  VFC  approach  to  control  robot  manipulators  for  contour  following  applications  in  [22] ,  and 
more  recently,  Dee  and  Li  used  VFC  to  achieve  passive  bilateral  teleoperation  of  robot  manipulators 
in  [18].  The  authors  of  [20]  utilized  a  passive  VFC  approach  to  develop  a  force  controller  for  robot 
manipulator  contour  following  applications.  Other  relevant  work  utilizing  VFC  approaches  are 
given  in  the  following:  [21]  and  [23].  Yamakita  et  al.  investigated  the  application  of  passive  VFC 
to  cooperative  mobile  robots  and  cooperative  robot  manipulators  in  [33]  and  [34],  respectively. 
Typically,  VFC  is  based  on  a  nonlinear  control  approach  where  exact  model  knowledge  of  the 
system  dynamics  are  required.  Motivated  by  the  desire  to  account  for  uncertainty  in  the  robot 
dynamics,  Cervantes  et  al.  developed  a  robust  VFC  in  [4],  Specifically,  in  [4]  a  proportional- 
integral  controller  was  developed  that  achieved  semiglobal  practical  stabilization  of  the  velocity 
field  tracking  errors  despite  uncertainty  in  the  robot  dynamics.  From  a  review  of  VFC  literature,  it 
can  also  be  determined  that  previous  research  efforts  have  focused  on  ensuring  the  robot  tracks  the 
velocity  field,  but  no  development  has  been  provided  to  ensure  the  link  position  remains  bounded. 
The  result  in  [4]  acknowledged  the  issue  of  boundedness  of  the  robot  position;  however,  the  issue  is 
simply  addressed  by  an  assumption  that  the  following  norm 


q(  0)  + 


d(q(a))da 


(1) 


yields  globally  bounded  trajectories,  where  q(t )  denotes  the  position,  and  0 ( ■ )  denotes  the  velocity 
field. 

In  addition  to  VFC,  some  task  objectives  are  motivated  by  the  need  to  follow  a  trajectory  to 
a  desired  goal  configuration  while  avoiding  known  obstacles  in  the  configuration  space.  For  this 
class  of  problems,  it  is  more  important  for  the  robot  to  follow  an  obstacle  free  path  to  the  desired 
goal  point  than  it  is  to  meet  a  time-based  requirement.  Numerous  researchers  have  investigated 
algorithms  to  address  this  motion  control  problem.  A  comprehensive  summary  of  techniques  that 
address  the  classic  geometric  problem  of  constructing  a  collision-free  path  and  traditional  path 
planning  algorithms  is  provided  in  Section  9,  “Literature  Landmarks”,  of  Chapter  1  of  [16].  Since 
the  pioneering  work  by  Khatib  in  [11],  it  is  clear  that  the  construction  and  use  of  potential  functions 
has  continued  to  be  one  of  the  mainstream  approaches  to  robotic  task  execution  among  known 
obstacles.  In  short,  potential  functions  produce  a  repulsive  potential  field  around  the  boundary  of 
the  robot  task-space  and  obstacles  and  an  attractive  potential  field  at  the  goal  configuration.  A 
comprehensive  overview  of  research  directed  at  potential  functions  is  provided  in  [16].  One  criticism 
of  the  potential  function  approach  is  that  local  minima  can  occur  that  can  cause  the  robot  to  “get 
stuck”  without  reaching  the  goal  position.  Several  researchers  have  proposed  approaches  to  address 
the  local  minima  issue  (e.g.,  see  [1],  [2],  [5],  [12],  [32]).  One  approach  to  address  the  local  minima 
issue  was  provided  by  Koditschek  in  [13]  for  holonomic  systems  (see  also  [14]  and  [27])  that  is  based 
on  a  special  kind  of  potential  function,  coined  a  navigation  function,  that  has  a  refined  mathematical 
structure  which  guarantees  a  unique  minimum  exists.  By  leveraging  from  previous  results  directed 
at  classic  (holonomic)  systems,  more  recent  research  has  focused  on  the  development  of  potential 
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function-based  approaches  for  nonholonomic  systems.  For  a  review  of  this  literature  see  [3],  [7],  [8], 
[9],  [15],  [17],  [23],  [25],  [27],  [30],  and  [31], 

The  aim  of  this  paper  is  to  illustrate  how  an  example  adaptive  controller  (e.g.,  the  benchmark 
adaptive  tracking  controller  presented  in  [28])  can  be  modified  to  incorporate  trajectory  planning 
techniques  with  the  controller.  To  this  end,  two  adaptive  controllers  are  developed.  The  first 
controller  focuses  on  the  VFC  problem.  Specifically,  the  benchmark  adaptive  controller  given  in 
[28]  is  modified  to  yield  VFC  in  the  presence  of  parametric  uncertainty.  The  contribution  of  the 
development  is  that  velocity  held  tracking  is  achieved  by  incorporating  a  norm  squared  gradient 
term  in  the  control  design  that  is  used  to  prove  the  link  positions  are  bounded  through  a  Lyapunov- 
analysis  rather  than  by  an  assumption.  In  lieu  of  the  assumption  in  (1),  the  VFC  development 
is  based  on  the  selection  of  a  velocity  held  that  is  first  order  differentiable,  and  that  a  first  order 
differentiable,  nonnegative  function  V(q)  G  R  exists  such  that  the  following  inequality  holds 

^t%)<-73(ll4ll)+<o  (2) 

where  9Vq^  denotes  the  partial  derivative  of  V(q )  with  respect  to  q(t),  73(-)  6  R  is  a  class  /C 
function2,  and  (0  G  R  is  a  nonnegative  constant.  That  is,  in  lieu  of  the  assumption  in  (1)  this 
paper  introduces  a  stability-based  condition  on  the  velocity  held.  It  is  interesting  to  note  that 
the  velocity  held  described  in  the  experimental  results  provided  in  [4]  can  be  shown  to  satisfy  the 
stability-based  condition  in  (2)  (see  Appendix  A  for  proof).  As  an  extension  to  the  VFC  problem, 
a  navigation  function  is  incorporated  with  the  benchmark  adaptive  controller  in  [28]  to  track  a 
reference  trajectory  that  yields  a  collision  free  path  to  a  constant  goal  point  in  an  obstacle  cluttered 
environment  with  known  obstacles. 

This  paper  is  organized  as  follows.  In  Section  2,  the  dynamic  model  for  a  robot  manipulator  is 
provided.  In  Section  3,  the  VFC  development  is  presented,  including  a  two-part  stability  analysis. 
The  first  analysis  proves  that  if  a  velocity  held  tracking  signal  is  square  integrable  then  the  link 
position  is  globally  uniformly  bounded  (GUB).  The  second  analysis  proves  that  the  velocity  held 
tracking  signal  is  square  integrable,  all  the  system  states  are  bounded,  and  that  the  velocity  held 
tracking  error  converges  to  zero  despite  parametric  uncertainty  in  the  dynamic  model.  Experimental 
results  based  on  the  velocity  held  presented  in  [4]  are  provided  to  demonstrate  proof  of  validation 
of  the  VFC  approach.  In  Section  4,  a  navigation  function  based  trajectory  planning  and  control 
development  is  presented,  along  with  the  stability  analysis.  This  analysis  proves  that  a  backstepping 
signal  is  square  integrable,  all  the  system  states  are  bounded,  and  that  the  robot  manipulator  will 
track  an  obstacle  free  path  to  a  goal  point,  despite  parametric  uncertainty  in  the  dynamic  model. 
Experimental  results  for  the  adaptive  navigation  function  controller  is  provided  to  demonstrate 
proof  of  validation  of  the  approach.  Concluding  remarks  are  provided  in  Section  5. 

2  System  Model 

The  mathematical  model  for  an  n-DOF  robotic  manipulator  is  assumed  to  have  the  following  form 

M(q)q  +  Vm(q,q)q  +  G(q)  =  t.  (3) 

In  (3),  q(t),  q(t),  q(t )  G  Mn  denote  the  link  position,  velocity,  and  acceleration,  respectively,  M(q)  G 
®mxn  represents  the  positive-dehnite,  symmetric  inertia  matrix,  Vm(q,q)  G  Mnxn  represents  the 

2  A  continuous  function  a  :  [0,  a)  — >  [0,  oo)  is  said  to  belong  to  class  K  if  it  is  strictly  increasing  and  a(0)  =  0  [10]. 
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centripetal-Coriolis  terms,  G(q)  G  BT  represents  the  known  gravitational  vector,  and  r(t)  G  Mn 
represents  the  torque  input  vector.  The  system  states,  q{t )  and  q{t)  are  assumed  to  be  measurable.  It 
is  also  assumed  that  M  (g),  Vm  (g,  g),  and  G  (g)  G  Coo  provided  g(f),  g(f)  G  Coo ■  The  dynamic  model 
in  (3),  exhibits  the  following  properties  that  are  utilized  in  the  subsequent  control  development  and 
stability  analysis. 

Property  1:  The  inertia  matrix  can  be  upper  and  lower  bounded  by  the  following  inequalities  [19] 

mi  ||5||2  <  ?M(q)f,  <  m2(q)  ||S||2  V£  e  R”  (4) 

where  mi  is  a  positive  constant,  m2(-)  is  a  positive  function,  and  ||-||  denotes  the  Euclidean 
norm. 

Property  2:  The  inertia  and  the  centripetal-Coriolis  matrices  satisfy  the  following  relationship 
[19] 

f  Qitf  (?)  -  vm(q,  q)  V  =  0  V(  €  K"  (5) 

where  M (q)  represents  the  time  derivative  of  the  inertia  matrix. 

Property  3:  The  robot  dynamics  given  in  (3)  can  be  linearly  parameterized  as  follows  [19] 

Y (Q,  Q,  Q)6  -  M(q)q  +  Vm(q,  q)q  +  G(q)  (6) 

where  6  G  Mp  contains  constant  system  parameters,  and  Y (g,  g,  g)  G  Mnxp  denotes  a  regression 
matrix  composed  of  q(t),  q(t),  and  q(t). 


3  Adaptive  VFC 

3.1  Control  Objective 

As  described  previously,  many  robotic  tasks  can  be  effectively  encapsulated  as  a  velocity  field.  That- 
is,  the  velocity  field  control  objective  can  be  described  as  commanding  the  robot  manipulator  to 
track  a  velocity  field  that  is  defined  as  a  function  of  the  current  link  position.  To  quantify  this 
objective,  a  velocity  field  tracking  error,  denoted  by  rj^t)  G  Rn,  is  defined  as  follows 

Vi(t)  =  q(t)  ~  #(q)  (7) 

where  $(•)  G  IT  denotes  the  velocity  field.  To  achieve  the  control  objective,  the  subsequent  devel¬ 
opment  is  based  on  the  assumption  that  q(t)  and  q(t)  are  measurable,  and  that  0(q)  and  its  partial 
derivative  e  HT,  are  assumed  to  be  bounded  provided  q(t)  G  Coo- 

3.2  Benchmark  Control  Modification 

To  develop  the  open- loop  error  dynamics  for  rj^t),  we  take  the  time  derivative  of  (7)  and  premultiply 
the  resulting  expression  by  the  inertia  matrix  as  follows 

M{q)f] i  =  —Vm(q,  q)q  —  G(q)  +  r  +  Vm(q,  g)'d(g)  (8) 

~Vm(q,  g)^(g)  -  M(g)^^-g 
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where  (3)  was  utilized.  From  (7),  the  expression  in  (8)  can  be  rewritten  as  follows 


=  ~Vm(q,  q)v i  -  Yi(q,  q)Q  +  r  (9) 

where  9  was  introduced  in  (6)  and  Y\  (q,q)  G  Mnxp  denotes  a  measurable  regression  matrix  that  is 
defined  as  follows 


Yi(q,  q)0  4  M(q)^q  +  Vm(q,  q)0(q)  +  G(q).  (10) 

Based  on  the  open-loop  error  system  in  (9),  a  number  of  control  designs  could  be  utilized  to  ensure 
velocity  held  tracking  (i.e. ,  ||r/1(f)||  — »  0)  given  the  assumption  in  (1).  Motivated  by  the  desire 
to  eliminate  the  assumption  in  (1),  a  norm  squared  gradient  term  is  incorporated  in  an  adaptive 
controller  introduced  in  [28]  as  follows 


T(t)± 


rj1  +  Y1(q,q)91 


(11) 


where  K  G  Wnxn  is  a  constant,  positive  definite  diagonal  matrix,  In  G  Mnxn  is  the  standard  n  x  n 
identity  matrix,  and  dVg ^  was  introduced  in  (2).  In  (11),  9  fit)  G  Mp  denotes  a  parameter  estimate 
that  is  generated  by  the  following  gradient  update  law 


91(t)  =  -T1Y1T(q,q)m  (12) 

where  Fi  G  M.pxp  is  a  constant,  positive  definite  diagonal  matrix.  After  substituting  (11)  into  (9), 
the  following  closed-loop  error  system  can  be  obtained 


M(q)rj  i 


-Vm(q,  q)Vi 


Yi(q,q)6 1  - 


dV(q ) 
dq 


where  the  parameter  estimation  error  signal  9\{t)  G  M.p  is  dehned  as  follows 

W)  =  9 -h. 


(13) 


(14) 


Remark  1  It  is  required  for  the  selection  of  a  particular  d(q)  and  V(q),  that  the  inequality  as 
defined  in  (2)  must  hold.  In  the  event  that  this  condition  does  not  hold,  the  tracking  objective  is  not 
guaranteed  as  described  by  the  subsequent  stability  analysis. 


Remark  2  While  the  control  development  is  based  on  a  modification  of  the  adaptive  controller 
introduced  in  [28],  the  norm,  squared  gradient  term  could  also  be  incorporated  in  other  benchmark 
controllers  to  yield  similar  results  ( e.g .,  sliding  mode  controllers) . 


3.3  Stability  Analysis 

To  facilitate  the  subsequent  stability  analysis,  the  following  preliminary  theorem  is  utilized. 
Theorem  1  Let  V(t)  Gi  denote  the  following  nonnegative,  continuous  differentiable  function 

V(t)  ^  V{q)  +  P(t) 
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where  V(q)  G  M.  denotes  a  nonnegative,  continuous  differentiable  function  that  satisfies  (2)  and  the 
following  inequalities 


0  <  Tidkll)  <  V(q)  <  Tadkll) 

where  71(-),  72(-)  are  class  /C  functions,  and  P(t )  G  M  denotes  the  following  nonnegative,  continuous 
differentiable  function 


P(t)  =  7  —  /  £2(cr)da 
Jt0 


(15) 


where  7  G  1  is  a  positive  constant,  and  s(t)  G  M  is  defined  as  follows 

|  9V  ( q ) 


A 

£  = 


dq 


(16) 


If  sit)  is  a  square  integrable  function,  where 


s2(a)da  <  7, 


I  to 


and  if  after  utilizing  (1),  the  time  derivative  ofV(t)  satisfies  the  following  inequality 


^W<-7s(lkll)+eo  (17) 

where  73(g)  is  the  class  /C  function  introduced  in  (2),  and  (0Gi  denotes  a  positive  constant,  then 
q(t )  is  global  uniformly  bounded. 

Proof:  The  time  derivative  of  V (t)  can  be  expressed  as  follows 


a  .  .  dViq )  .  .  dViq ) 

V  (; t )  = 


<9g 


<9g 


where  (7)  and  (15)  were  utilized.  By  exploiting  the  inequality  introduced  in  (2)  and  the  definition 
for  s{t)  provided  in  (16),  the  following  inequality  can  be  obtained 


v  (t)  <  -73(IMI)  +  Co  +  [£(t)  -  £2(t)]  ■  (!8) 

After  completing  the  squares  on  the  bracketed  terms  in  (18),  the  inequality  introduced  in  (17)  is 
obtained  where 

Co  —  Co  +  |- 

Hence,  if  s(t)  G  then  Lemma  1  in  Appendix  B  can  be  used  to  prove  that  q{t)  is  GUB.  ■ 

In  the  following  analysis,  we  first  prove  that  e(t)  G  C2-  Based  on  the  conclusion  that  e(t)  G  C2, 
the  result  from  Theorem  1  is  utilized  to  ensure  that  q{t)  is  bounded  under  the  proposed  adaptive 
controller  given  in  (11)  and  (12). 


Theorem  2  The  adaptive  VFC  given  in  (11)  and  (12)  yields  global  velocity  field  tracking  in  the 
sense  that 


ll^iWII  "^0  as  t  ^  00. 


(19) 
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Proof:  Let  V\  (t)  G  M  denote  the  following  nonnegative  function 


Vi  =  ln[Mm  +  iflhr'Si. 

After  taking  the  time  derivative  of  (20)  the  following  expression  can  be  obtained 


(20) 


V)  -  -ri‘,  Y1(g,q)'e1+  K  + 


dV(q ) 


dq 


ihMirr1 0i 


(21) 


where  (5)  and  (13)  were  utilized.  After  utilizing  the  parameter  update  law  given  in  (12),  the 
expression  given  in  (21)  can  be  rewritten  as  follows 


(22) 


The  expressions  given  in  (16),  (20),  and  (22)  can  be  used  to  conclude  that  rh(t),  9\{t)  G  Coo  and 
rh(t),£(t)  G  £2-  Based  on  the  fact  that  e(t)  G  £2,  the  results  from  Theorem  1  can  be  used  to  prove 
that  q(t)  G  Since  q(t)  G  Coo,  the  assumption  that  i9(q)  and^p  G  Coo  can  be  used  to  conclude 

that  q  (t)  G  £oo,  where  the  expression  in  (7)  was  utilized.  Based  on  the  fact  that  9i(t)  G  £oo,  the 
expression  in  (14)  can  be  used  to  prove  that  9i(t)  G  Coo ■  Based  on  the  facts  that  M  (q),  Vm  (q,  q), 
G  ( q )  G  £oo  along  with  the  facts  that  q  (t) ,  q  (' t ),  i9(q),  G  Coo,  and  (10)  can  be  used  to  prove 

that  y\(q,  q)  G  £ 00 ■  The  facts  that  Y\(q,  q),  9i(t),  rj^t)  G  Coo  and  the  fact  that  dv^  G  Coo  (based  on 
Theorem  1)  can  be  used  along  with  (11)  to  prove  that  r(t)  G  Coo ■  Based  on  the  previous  bounding 
statements,  the  expression  given  in  (13)  can  be  used  to  prove  that  rji(t)  G  Coo-  Given  that  rh(t), 
rh(t)  G  Coo  and  rji(t)  G  £2,  Barbalat’s  Lemma  [28]  can  now  be  utilized  to  prove  (19).  ■ 


4  Navigation  Function  Control  Extension 

4.1  Control  Objective 

The  objective  in  this  extension  is  to  navigate  a  robot’s  end-effector  along  a  collision-free  path  to  a 
constant  goal  point,  denoted  by  q*  G  V,  where  the  set  D  denotes  a  free  configuration  space  that  is 
a  subset  of  the  whole  configuration  space  with  all  configurations  removed  that  involve  a  collision 
with  an  obstacle,  and  q*  G  M"‘  denotes  the  constant  goal  point  in  the  interior  of  T>.  Mathematically, 
the  primary  control  objective  can  be  stated  as  the  desire  to  ensure  that 

q(t)  — ►  q*  as  t  — >  00  (23) 

where  the  secondary  control  is  to  ensure  that  q(t)  G  T>.  To  achieve  these  two  control  objectives,  we 
define  (p  (q)  G  I  as  a  function  (p(q)  :  T>  — >[0,  1]  that  is  assumed  to  satisfy  the  following  properties: 

•  PI)  The  function  p  ( q )  is  first  order  and  second  order  differentiable  Morse  function  [14]  (i.e. , 
-§-qv{q) and  -§-q  (jfc<p(q))  exist  on  V). 

•  P2)  The  function  p  ( q )  obtains  its  maximum  value  on  the  boundary  of  T>. 

•  P3)  The  function  (q)  has  a  unique  global  minimum  at  q  (t)  =  q* . 
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P4)  If  -§^(p{q)  =  0  then  q  (t)  =  q*. 


Based  on  (23)  and  the  above  definition,  an  auxiliary  tracking  error  signal,  denoted  by  r)2(t)  £  Mn, 
can  be  defined  as  follows  to  quantify  the  control  objective 

Vi(t)  =  q(t)  +  S7<p(q)  (24) 

where  \/tp  (q)  =  -j^cp  (q)  denotes  the  gradient  vector  of  ip  ( q )  defined  as  follows 


v<p(q)  = 


dip  dip 

dqi  dq2 


dip 

dqn 


l  T 


(25) 


Remark  3  As  discussed  in  [27],  the  construction  of  the  function  ip  (q),  coined  a  navigation  function, 
that  satisfies  all  of  the  above  properties  for  a  general  obstacle  avoidance  problem .  is  nontrivial.  Indeed, 
for  a  typical  obstacle  avoidance,  it  does  not  seem,  possible  to  construct  ip  (q)  such  that  ip  (q)  =  0 
only  at  q[t )  =  q*.  That  is,  as  discussed  in  [27],  the  appearance  of  interior  saddle  points  (i.e., 
unstable  equilibria)  seems  to  be  unavoidable;  however,  these  unstable  equilibria  may  have  minimal 
impact  in  practice.  That  is,  ip  (q)  can  be  constructed  as  shown  in  [27]  such  that  only  a  “few”  initial 
conditions  will  result  in  convergence  to  the  unstable  equilibria. 


Remark  4  The  two  control  developments  presented  in  the  paper  appear  to  be  mathematically  similar 
(i.e.  (7)  and  (24)),  but.  the  control  objectives  are  very  different..  The  VFC  objective  is  to  achieve 
robot,  end-effector  velocity  tracking  with  a  desired  trajectory  generated  by  a  velocity  field,  id  (q),  hence, 
there  is  no  goal  point. .  The  navigation  function  control  development,  utilizes  a  special  function  (p  (q) , 
that  has  specific  properties  such  that,  the  robot’s  end-effector  finds  a  collision  free  path  to  a  known 
goal  point.,  q*  and  stops.  Each  signal,  id  (q)  for  the  VFC  development,  and  <p  (q)  for  the  navigation 
function  control  development,  must,  meet  a  set.  of  qualifying  conditions  (i.e.  the  inequality  of  (2)  for 
the  VFC  and  PI  -  P4  for  the  navigation  function  control),  but.  these  conditions  are  not.  the  sam.e, 
therefore  the  two  objectives  are  very  different. 


4.2  Benchmark  Control  Modification 

To  develop  the  open- loop  error  dynamics  for  rj2(t),  we  take  the  time  derivative  of  (24)  and  premul¬ 
tiply  the  resulting  expression  by  the  inertia  matrix  as  follows 

Mi]2  =  -Vm(q,  q)v 2  +  Y2(q,  q)0  +  r.  (26) 

where  (3)  and  (24)  were  utilized.  In  (26),  the  linear  parameterization  Y2(q,q)6  is  defined  as  follows 

Y2{q,  q)6  4  M(q)f(q,  q)  +  Vm(q,  q)  V  <p(q)  ~  G{q )  (27) 

where  Y2(q,q)  £  Mnxm  denotes  a  measurable  regression  matrix,  6  £  Rm  was  introduced  in  (6),  and 
the  auxiliary  signal  f(q,  q)  £  Rn  is  defined  as 

/M)  =  ^(w  (q)) 

=  H(q)q 
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(28) 


where  the  Hessian  matrix  H(q )  G  Mnxn  is  defined  as  follows 


H{q)  = 


'  d2ip 

d2(p 

d2cp 

dq\ 

<3>V 

dqxdq2 

d2tp 

dqidqn 

d2p 

dq2dq1 

dq2 

dq2dqn 

d2ip 

d2p 

_  dqndq1 

9q2n 

Based  on  (26)  and  the  subsequent  stability  analysis,  the  following  adaptive  controller  introduced  in 
[28]  can  be  utilized 

t  =  -krj2-Y2(q,q)92  (29) 


where  lei  is  a  positive  constant  gain,  and  92{t )  G  Rp  denotes  a  parameter  update  law  that  is 
generated  from  the  following  expression 


92(t)±r2Y?(q,q)ri2  (30) 

where  F2  G  Mmxm  is  a  positive  definite,  diagonal  gain  matrix.  Note  that  the  trajectory  planning 
is  incorporated  in  the  controller  through  the  gradient  terms  included  in  (27)  and  (28).  After 
substituting  (29)  into  (26)  the  following  closed  loop  error  systems  can  be  obtained 

Mri2  =  -Vm(q,  q)g2  -  kr]2  +  Y2(q,  q)92  (31) 

where  92{t )  G  Mp  is  defined  as  follows 

~e2{t)^  9  -  92.  (32) 


4.3  Stability  Analysis 


Theorem  3  The  adaptive  controller  given  in  (29)  and  (30)  ensures  that  the  robot  manipulator 
tracks  an  obstacle  free  path  to  the  unique  goal  configuration  in  sense  that 


q(t )  — >  q*  as  t  — >  oo 

provided  the  control  gain  k  introduced  in  (29)  is  selected  sufficiently  large. 
Proof:  Let  V2(q,  q2,92)  G  M  denote  the  following  nonnegative  function 


(33) 


where  7  G  1  is  an  adjustable,  positive  constant.  After  taking  the  time  derivative  of  (33)  the 
following  expression  can  be  obtained 


v2  =  [VT(q)Y  Q  +  Ihl  +  *2 (q,  #2)  -  1~9T2. r2  1  92 
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where  (5),  (25),  (31),  and  (32)  were  utilized.  By  utilizing  (24),  (30),  the  following  expression  can 
be  obtained 

v2  =  -  ||w(?)ll2  -  jk  ||%||2  +  lvc(q)fv2- 

The  expression  above  can  be  further  simplified  as  follows 

V2  <  —  II  V^(?)||2  -  (7*  -  2)  IM2  (34) 

where  the  following  upper  bound  was  utilized 

[w(?)]T»7: 2  <  \  llv^WII2  +  2||%H2- 

Provided  k  is  selected  sufficiently  large  to  satisfy 

k  >  — ,  (35) 

7 

it  is  clear  from  (4),  (33),  and  (34)  that 

0  <  <p(q(t))  +  7C (q,  t )  <  <p(g(0))  +  7C(?(0),  0)  (36) 

where  ( (q,t)  G  M  is  dehned  as 

((«,*)=  [^H>)2WII2  +  Wr2-1}||e2W||2].  (37) 

From  (32),  (36),  and  (37)  it  is  clear  that  r)2(t),  <p(q),  O2 (t),  02 (t)  E  Coo ■  Let  the  region  Vq  be  dehned 
as  follows 

^0  =  {q{t) |0  <  <p{q(t))  <  <p{q( 0))  +  7C(?(0),  0)}  .  (38) 

Hence;  (33),  (34)  and  (36)  can  be  utilized  to  show  that  q(t)  E  Vq  provided  q(0)  E  V 0  (i.e. ,  q(t)  E  V0 
Vg( 0)  E  T>0).  Based  on  property  PI  given  above,  we  now  know  that  \Jtp(q)  £  C^  Vg( 0)  E  V 0. 
Since  rj2 E  Coo  Vg( 0)  E  V 0,  (24)  can  be  used  to  conclude  that  q(t)  E  Coo  Vg( 0)  E  V0-, 
hence,  property  PI)  and  (28)  can  be  used  to  conclude  that  f(q,q )  G  Coo  Vg( 0)  G  V0.  Based  on  the 
fact  that  M  (q),  Vm(q,q),  G  (q)  E  Coo  Vg( 0)  G  T>0  along  with  the  facts  that  S/q?(q),  f(q,q )  G  Coo 
Vg( 0)  G  Vq  can  be  used  along  with  (27)  to  prove  that  Y2(q,q)  E  Coo  Vg( 0)  G  T>0.  Based  on 
the  facts  that  r)2(t),  Y2(q,q),  02(t)  E  Coo  Vg( 0)  G  V0  can  be  used  along  with  (29)  to  prove  that 
r(t)  E  Coo  Vg( 0)  G  T>0.  Based  on  the  previous  boundedness  statements,  (31)  can  be  used  to  show 
that  f] 2{t)  E  £ooVg(0)  G  Vq,  and  hence,  S/q?(q),  rj2( t )  are  uniformly  continuous  Vg(0)  G  D0.  From 
(34),  it  can  also  be  determined  that  S/q?(q),  T}2(t)  £  C2  Vq(0)  E  V0.  From  these  facts,  Barbalat’s 
Lemma  [28]  can  be  used  to  show  that  Vc(q)>  V2 (t)  0  as  t  — >  00  Vg( 0)  G  V0.  Since  S/q?(q)  — >  0 

,  property  P4)  given  above  can  be  used  to  prove  that  q{t)  — >  q*  as  t  — >•  00  Vg( 0)  G  Vq.  To  ensure 
that  q{t)  will  remain  in  a  collision-free  region,  we  must  account  for  the  effects  of  the  7C(7,(0),0) 
term  introduced  in  the  definition  of  the  region  Vq  given  in  (38).  To  this  end,  we  first  define  the 
region  Vi  as  follows 

£>i  -  {q(t)\ 0  <  <p(q(t))  <  1}  (39) 
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where  T> 4  denotes  the  largest  collision-free  region,  which  is  based  on  the  dehnition  of  the  function 
<p  (q)  :  V  — >[0,  1].  It  is  now  clear  from  (38)  and  (39)  that  if  the  weighting  constant  7  is  selected 
sufficiently  small  to  satisfy 


¥»(?(0))  +  7C(?(0),0)<l,  (40) 

hence  making  the  upper  bound  of  V\  greater  than  the  upper  bound  of  V0,  then  'D{]  C  V 4,  and 
hence,  the  robot  manipulator  tracks  an  obstacle  free  path.  ■ 


5  Experimental  Verification 

Experimental  results  were  obtained  by  implementing  the  adaptive  VFC  and  the  navigation  function 
controller  on  a  Barrett  Whole  Arm  Manipulator  (WAM).  The  experimental  testbed  and  results  from 
implementing  the  controllers  are  provided  in  the  following  sections. 

5.1  Experimental  Setup 

The  WAM  testbed  depicted  in  Fig.  1  was  utilized  to  implement  the  VFC  and  the  navigation 
function  controller.  For  simplicity,  5  links  of  the  robot  were  locked  at  a  fixed,  specified  angle  during 
the  experiment,  and  the  remaining  links  of  the  manipulator  were  used  to  enable  the  manipulator  to 
move  along  a  planar  trajectory.  Specifically,  a  joint-space  proportional  derivative  (PD)  controller 
was  utilized  to  servo  the  WMR  to  the  following  initial  joint  configuration  for  the  adaptive  VFC 
experiment  (in  [deg]) 


g(0)=[  0  90  -90  60  90  20  0  ]T 

and  to  the  following  joint  configuration  for  the  navigation  function  experiment  (in  [deg]) 

g(0)=[  -58.84  90  90  140.72  11.5  84.5  0  ]T. 


Once  the  WAM  was  servoed  to  the  initial  joint  configuration,  links  2,  3,  5,  6  and  7  were  locked, 
resulting  in  a  planar  configuration  with  links  1  and  4  (see  Fig.  1).  The  resulting  forward  kinematics 
and  manipulator  Jacobian  for  the  planar- WAM  are  given  as  follows 


Xi 

£1  cos(gi)  +  £4  cos(gi  +  q4) 

-  x2  . 

£4  sin(qi)  +  £4  sin(g1  +  g4) 

J{q) 


—l\  sin(gi)  -  £4  sin(gi  +  q4)  —l4  sin  (<71  +  q4) 

£ 1  cos(gi)  +  £4  cos(gi  +  q4)  £4  cos(gi  +  q4) 


(41) 


where  £\  =  0.558  [m]  and  £4  =  0.291  [m].  The  dynamics  of  the  planar- WAM  can  be  expressed  in 
the  following  form  [29] 


r  = 


+ 


Mu  M12 

qi 

M2i  M22 

..  ^  .. 

'4,  0 

Qi 

.  °  fd4  _ 

.  94  . 

+ 


Knn 

Pm2i 


147112 

14n22 


(42) 
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In  (42),  the  elements  of  the  inertia  and  centripetal-Coriolis  matrices  are  given  as 


Mn  =  Pi  +  2  p2cos(q4) 
Mi2  =  p3  +p2cos(q4) 
M2 1  =  p3  +P2Cos(g4) 
M22  =  p3 


Knn  =  -'P'lSinipuVu 

vmi2  =  -P2sin(q4)q\  -  P2sin(q4)q4 

Vm21  =  p2sin(q4)q1 

Vm22  =  0 

where  pi,  p2,  P:>,  denote  unknown  constant  inertial  parameters,  and  =  6.8  [Nm-s]  and  fd4  = 
3.8  [Nm-s].  The  gravitational  effects  are  not  included  in  (42)  due  to  the  plane  of  motion  of  the 
manipulator. 

The  links  of  the  WAM  manipulator  are  driven  by  brushless  motors  supplied  with  sinusoidal 
electronic  commutation.  Each  axis  has  encoders  located  at  the  motor  shaft  for  link  position  mea¬ 
surements.  Since  no  tachometers  are  present  for  velocity  measurements,  link  velocity  signals  are 
calculated  via  a  filtered  backwards  difference  algorithm.  An  AMD  Athlon  1.2GHz  PC  operating 
QNX  6.2.1  RTP  (Real  Time  Platform),  a  real-time  micro-kernel  based  operating  system,  hosts  the 
control,  detection  and  identification  algorithms  which  were  written  in  “C-|— 1-”.  Qmotor  3.0  [24],  was 
used  to  facilitate  real  time  graphing,  data  logging  and  on-line  gain  adjustment.  Data  acquisition 
and  control  implementation  were  performed  at  a  frequency  of  1.0  [kHz]  using  the  ServoToGo  I/O 
board. 


Figure  1:  Front  view  of  the  experimental  setup 
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5.2  Experimental  Results 

5.2.1  Adaptive  VFC  Experiment 

The  following  task-space  velocity  field  for  a  planar,  circular  contour  was  utilized  for  the  experiment 

[4] 


$(x) 


-2  K(x)f(x) 


(xi  -  Xcl) 

{x2  ~  Xc2) 


+  2c(x) 


~{X2  ~  Xc2) 

-  Xci) 


(43) 


where  xci  =  0.54  [m]  and  xc2  =  0.04  [m] 
c(x)  G  M  are  defined  as  follows  [4] 

f(x)  = 
K(x)  = 


c(x)  = 


denote  the  circle  center,  and  the  functions  f(x),  K(x),  and 


(xi  -  xci)2  +  (x2  -  xc2)2  -  r 
h* 


sTs W) 


df(x) 

dx 


+  6 

Co  exp  f  -/'  v'/-(o  j 


9f(x) 

dx 

(44) 


In  (44),  rQ  =  0.2  [m]  denotes  the  circle  diameter,  and  the  parameters  e  =  0.005  [m3],  //  =  20  [m_1], 
k*  =  0.25  [ms-1],  and  c0  =  0.25  [ms-1]  were  selected  according  to  [4],  The  task-space  velocity  field 
is  depicted  in  Fig.  2.  The  development  in  Appendix  A  indicates  that  the  velocity  field  in  (43) 
satisfies  the  condition  given  in  (2).  To  implement  the  adaptive  VFC  given  in  (11)  and  (12),  the 
task-space  velocity  field  is  transformed  into  a  joint-space  velocity  field  as  follows  id(q)  =  J~1(q),d(x). 

The  following  values  were  recorded  as  follows 


K  =  diag( 25, 15)  T  =  diag( 3, 1,  5) 

where  diag(-)  denotes  a  diagonal  matrix  with  the  arguments  as  the  diagonal  entries,  the  resulting 
velocity  field  tracking  errors  are  given  in  Fig.  2.  Fig.  4  depicts  the  parameter  estimates,  and  Fig. 
5  depicts  the  control  torque  inputs. 


5.2.2  Adaptive  Navigation  Function  Control  Experiment 

For  the  navigation  function  control  experiment,  four  circular  obstacles  with  known  dimensions  were 
placed  in  known  locations  in  the  task-space  (see  Fig.  1).  The  actual  size  of  the  obstacles  and 
task-space  is  then  modified  in  the  algorithm  to  accommodate  for  the  term  7C(o,(0))0)  given  in 
(36)  and  (37)  (i.e. ,  the  configuration-space  was  reduced  to  ensure  obstacle  avoidance).  To  modify 
the  configuration-space  according  to  (36)  and  (37),  exact  knowledge  of  the  inertial  parameters  is 
required.  Since  these  parameters  are  unknown,  an  upper  bound  for  C(<?(0),  0)  was  utilized  based  on 
known  upper  bounds  for  the  inertial  parameters.  The  modifications  to  the  configuration-space  are 
depicted  in  Figure  6.  A  task-space  navigation  function  was  developed  to  encapsulate  the  obstacles 
and  the  task-space  boundary  as  follows 


(||x  -  x*||28  +  AA/^AA)1714 


(45) 
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Figure  2:  Desired  Trajectory 


Figure  3:  Velocity  field  tracking  errors 
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Figure  5:  Control  torque  inputs. 


where  x(t)  =  [xi(t),  x2(t)]T  G  M2  denote  the  actual  task-space  position  of  the  WAM  end-effector, 
x*  =  [x*,  X2]71  G  M2  denotes  the  task-space  goal  position.  In  (45),  the  boundary  function  /30(x )  G  M 
and  the  obstacle  functions  Pi(x),  P2(x),  /3s(x),  /34(x)  G  M  are  defined  as  follows 

Po  =  r%-  (xi  -  xiro)2  -  (x2  -  x2ro)2  (46) 

Pi  =  (xi  ~  xri)2  +  (x2  —  x2ri)2  —  r2 

p2  =  {xi  -xr2f  +  (x2  -X2r2)2  -r\ 

Pz  =  (xi  -  Xr3)2  +  (x2  -  x2rs)2  -  r\ 

Pa  =  {xi  ~  xr4)2  +  (x2  -  x2r4)2  -  r\. 

In  (46),  (xi  —  X\ri)  and  (x2  —  x2n)  where  i  —  0,1,  2,  3, 4  are  the  respective  centers  of  the  boundary 
and  obstacles,  and  xq,  r\:  r2,  r3,  r4  G  M  are  the  respective  radii  of  the  boundary  and  obstacles.  From 
(45)  and  (46)  it  is  clear  that  the  model-space  is  a  circle  that  excludes  four  smaller  circles  described 
by  the  obstacle  functions  Ppx),  P2(x),  P3(x),  /?4(x).  Based  on  the  known  location  and  size  of 
the  obstacles  and  task-space  boundary,  the  model-space  configuration  parameters  were  selected  as 
follows  (in  [m]) 

Xiro  =  0.5064  x2ro  =  —0.0275  r0  =  0.28 

xiri  =  0.63703  x2ri  =  0.11342  r4  =  0.03 

x\ r2  =  0.4011  x2r2  =  0.0735  r2  =  0.03 

X\ r3  =  0.3788  x2rs  =  —0.1529  r3  =  0.03 

X\ r4  =  0.6336  x2r4  =  —0.12689  r4  =  0.03. 

To  implement  the  navigation  function  based  controller  given  in  (29)  and  (30)  the  joint-space 
dynamic  model  given  in  (42)  was  transformed  to  the  task-space  as  follows  [6] 

T*  =  M*(x)  x  +  V^(x,x)x  +  F^x 

where 


where  J  1(q)  can  be  determined  from  (41)  as  follows 

cos  (gi  +  (u)  sin  (gi  +  g4) 

j- i(g\  =  h  sin  q4  t\  sin  q4 

i\  cos  q4  +  £4  cos  (qi  +  q4)  t\  sin  q4  +  £4  sin  (q4  +  g4) 

i4l4  sin  q4  £4£4  sin  q4 

After  adjusting  the  control  gains  to  ensure  the  gain  conditions  (35)  and  (40)  are  satisfied,  the 
following  values  were  recorded 


k  =  45  r2  =  diag  (0.02, 0.01,  0.01), 


the  resulting  actual  trajectory  of  the  WAM  end-effector  is  provided  in  Fig.  6.  Fig.  6  illustrates  that 
the  WAM  end-effector  avoids  the  actual  obstacles  as  it  moves  to  the  goal  point.  The  parameter 
estimates  and  control  torque  inputs  are  provided  in  Figs.  7  and  8,  respectively. 


16 


Figure  6:  Actual  trajectory  of  the  WAM  robot. 


Figure  7:  Parameter  estimates. 


Link  1 


Figure  8:  Control  torque  inputs. 


Remark  5  The  adaptive  control  results  achieved  in  Sections  3.3  and  4-3  proves  that  6 i  ( t )  and 

02  (t)  — > ►  0  as  t  — >  oo,  therefore  the  parameter  estimates  6  (t)  are  not  identified.  The  values  of  the 
estimates  that  are  reached  could  be  different  from,  one  experimental  run  to  another.  The  parameter 
estimates  may  not  become  constant  due  to  steady-state  tracking  errors. 

6  Conclusion 

Two  trajectory  planning  and  adaptive  tracking  controllers  are  presented.  The  benchmark  adaptive 
tracking  controller  by  Slotine  [28]  was  modified  to  achieve  velocity  field  tracking  in  the  presence 
of  parametric  uncertainty  in  the  robot  dynamics.  By  incorporating  a  norm  squared  gradient  term 
to  the  VFC,  the  boundedness  of  all  signals  can  be  proven  without  the  typical  assumption  that 
bounds  the  integral  of  the  velocity  field.  An  extension  was  then  provided  that  also  modifies  a 
standard  adaptive  controller  by  incorporating  a  gradient  based  term.  Using  standard  backstepping 
techniques,  a  Lyapunov  analysis  was  used  to  prove  that  a  navigation  function  could  be  incorporated 
in  the  control  design  to  ensure  the  robot  remained  on  an  obstacle  free  path  within  an  expanded 
configuration  space  to  reach  a  goal  configuration.  Experimentation  results  illustrated  proof  of 
validation  of  the  adaptive  VFC  and  navigation  function  tracking  controllers. 
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A  Experimental  Velocity  Field  Selection 

This  VFC  development  is  based  on  the  selection  of  a  velocity  held  that  is  first  order  differentiable, 
and  that  a  first  order  differentiable,  nonnegative  function  V(q)  G  M.  exists  such  that  the  following 
inequality  holds 

—  — 73(IMI)  +  Co  (47) 

where  9vd^  denotes  the  partial  derivative  of  V(q)  with  respect  to  q{t),  73(-)  6  R  is  a  class  /C 
function3,  and  (0  e  M  is  a  nonnegative  constant.  To  prove  that  the  velocity  held  in  (43)  and 
(44)  satishes  the  condition  in  (47),  let  V(x)  G  R.  denote  the  following  nonnegative,  continuous 
differentiable  function 

V(x)  =  ^ xTx .  (48) 

After  taking  the  time  derivative  of  (48)  and  substituting  (43)  for  x(t),  the  following  inequality  can 
be  developed 

V  =  xTtf(x)  <  -73(x)  +  Co  (49) 

where  73(a;)  and  Co  were  dehned  in  (47). 

To  prove  the  inequality  given  in  (49)  we  must  hnd  q3(x)  and  Co-  To  this  end,  we  rewrite  xTtt(x) 
as  follows 


xTd(x)  =  —2  K(x)f(x)xTx 


(50) 


where  (43)  has  been  utilized,  and  xc\  and  ay 2  of  (43)  and  (44)  are  set  to  zero  for  simplicity  and 
without  loss  of  generality.  By  substituting  (44)  into  (50)  for  K(x)  and  f(x),  the  following  expression 
can  be  obtained 


_ u*  I  ,T\  4  1  u*r 2  I  ,r|  2 

|  'E  |  rv q  I  Q  |  dj  | 

|(H2->7|H+C 


(51) 


After  utilizing  the  following  inequality 


x 


<«iwr  + 


1 

6 


3  A  continuous  function  a  :  [0,  a)  — *  [0, 00)  is  said  to  belong  to  class  K  if  it  is  strictly  increasing  and  ck(0)  =  0  [10]. 


21 


where  <5  G  R  is  a  positive  constant,  the  inequality  given  in  (49)  can  be  determined  from  (51)  where 


73(x)  = 


k*(l^r*6)\\x\\4 

(IMI2  +  ro)  IMI  +  ! 


and 


Co  — 


2  Krl 
6e  ' 


Provided  <5  is  selected  according  to  the  following  inequality 

6  <  4, 


then  73(a;)  can  be  shown  to  be  a  class  /C  function. 


B  Lemma  1 

Given  a  continuously  differentiable  function,  denoted  by  V (q).  that  satishes  the  following  inequalities 

0  <  Tidkll)  <  V(q)  <  72(lkll)  +Cfe  (52) 

with  a  time  derivative  that  satisfies  the  following  inequality 

V(q)  <  — 73(IMI)  +  Co,  (53) 

then  q{t )  is  GUB,  where  q^-),  72(');  7s(')  are  class  /C  functions,  and  £0,  G  R  denote  positive 
constants. 

Proof:1  Let  fi  G  M  be  a  positive  function  defined  as  follows 


G  —  73  1(Co)  >  0  (54) 

where  y3x  (•)  denotes  the  inverse  of  73(-)>  and  let  £>(0,f2)  denote  a  ball  centered  about  the  origin 
with  a  radius  of  fh  Consider  the  following  2  possible  cases. 

The  initial  condition  g(f0 )  lies  outside  the  ball  B( 0,  G)  as  follows 

f2  <  \\q(t0)\\  <  Gi  (55) 

where  hi  G  1  is  a  positive  constant.  To  facilitate  further  analysis,  we  define  the  operator  d(-)  as 
follows 


d(fii)  -  (7i 1  °  7a)(fii)  +  7i  X(C6)  >  0  (56) 

where  (7b1  0  72)  denotes  the  composition  of  the  inverse  of  (•)  with  72  (•)  (i.e.,  the  inverse  of  the 
function  r)l  (•)  is  applied  to  the  function  72  (•)).  After  substituting  the  constant  g^Gx)  into  q^-), 
the  following  inequalities  can  be  determined 

7i (d  (fii))  =  72(^i)  +  C&  >  72(lk(^o)||)  +  Cb  >  V  (. q(t0 ))  (57) 

4 This  proof  is  based  on  the  proof  for  Theorem  2.14  in  [26]. 
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where  the  inequalities  provided  in  (52)  and  (55)  were  utilized. 

Assume  that  q  (r)  €  R  for  t0  <  r  <  t  <  oo  lies  outside  the  ball  B( 0, 12)  as  follows 

12  <  || 5  (t)  || . 

From  (53)  and  (58),  the  following  inequality  can  be  determined 

V  ( q  (t))  <  -73 (12)  +  £0> 

and  hence,  from  the  definition  for  12  in  (54),  it  is  clear  that 

V(q(r))<  0. 


(58) 


(59) 


By  utilizing  (57)  and  the  result  in  (59),  the  following  inequalities  can  be  developed  for  some  constant 
At 


7i (d  (fii))  >  V  ( q(t0 ))  >  V  (q(r))  >  V  (q(r  +  At))  >  (||g(r  +  Ar)||)  .  (60) 

Since  7X(-)  is  a  class  /C  function,  (56)  and  (60)  can  be  used  to  develop  the  following  inequality 

h(t)\\  <  d(fii)  =  (7b1  °  72)(^i)  +  7b1  (^)  Vt  >  t0 

provided  the  assumption  in  (58)  is  satisfied.  If  the  assumption  in  (58)  is  not  satisfied,  then 

||g(b>||  <  12  =  7b1(^o)  V2>20-  (61) 


Hence,  q(t)  is  GUB  for  Case  A. 

The  initial  condition  q(to)  lies  inside  the  ball  B( 0, 12)  as  follows 

h(t0)\\  <  12  <  12i. 

If  q(t)  remains  in  the  ball,  then  the  inequality  developed  in  (61)  will  be  satisfied.  If  q(t)  leaves  the 
ball,  then  the  results  from  Case  A  can  be  applied.  Hence,  q(t)  is  GUB  for  Case  B.  ■ 
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